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ABSTRACT 

We present a high-precision J-band light curve for the Wolf-Rayet binary 
WR 20a, obtained as a sub-project of the Optical Gravitational Lensing Experi- 
ment. Rauw et al. have recently presented spectroscopy for this system, strongly 
suggesting extremely large minimum masses of 70.7 ± 4.0 M and 68.8 ± 3.8 M Q 
for the component stars of the system, with the exact values depending strongly 
on the period of the system. We detect deep eclipses of about 0.4 mag in the light 
curve of WR 20a, confirming and refining the suspected period of P = 3.686 days 
and deriving an inclination angle of i = 74.5 ± 2.0 deg. Using these photometric 
data and the radial velocity data of Rauw et al., we derive the masses for the 
two components of WR 20a to be 83.0 ± 5.0 M and 82.0 ± 5.0 M . Therefore, 
WR 20a is confirmed to consist of two extremely massive stars and to be the 
most massive binary known with an accurate mass determination. 

Subject headings: binaries: eclipsing - binaries: spectroscopic - stars: Wolf-Rayet 
- stars: fundamental parameters - stars: individual: WR 20a 



1 Based on observations obtained with the 1.3 m Warsaw telescope at Las Campanas Observatory, which 
is operated by the Carnegie Institute of Washington. 
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1. Introduction 

Measuring accurate masses for the most massive stars in our Galaxy and beyond is 
important for constraining star formation and stellar evolution theories, which have indirect 
implications for many objects that are not well understood, such as supernovae, gamma-ray 
bursts and Population III stars. The most massive candidates in the Milky Way are LBV 
1806-20 (Eikenberry et al. 2004) and the Pistol Star (Figer et al. 1998), which have inferred 
masses up to ~ 200 M Q . HDE 269810 in the LMC is another massive candidate ~ 150 M 
(Walborn et al. 2004). Finally, rj Carinae, one of the most studied massive stars has an 
estimated mass of ~ 120 M , yet it remains unknown whether it is a rapid rotator or a close 
binary system (see review by Davidson & Humphreys 1997). However, all these masses are 
only indirect estimates and thus have large uncertainties associated with them. The only 
direct way of measuring accurate masses of distant stars is through double lined spectroscopic 
binary systems with eclipses present in their light curves. The limited number of existing 
mass measurements for massive stars (> 20 M ) can be explained by their intrinsic rarity 
and the observationally demanding process of discovering them in eclipsing binary systems. 

Until recently, the most massive stars ever weighed in binaries were: R136-38 (03V+06V) 
in the LMC with a primary mass of 56.9±0.6 M (Massey et al. 2002), WR 22 (WN7+abs + O), 
with a minimum primary mass of 55.3 ±7.3 M (Rauw et al. 1996; Schweickhardt et al. 1999) 
and Plaskett's star with a minimum primary mass of 51 M (Bagnuolo et al. 1992). How- 
ever, WR 20a, a Wolf-Rayet (WR) binary in the compact cluster Westerlund 2, seems to be 
the new record holder. Rauw et al. (2004) obtained spectroscopy for WR 20a and measured 
extremely large minimum masses of 70.7 ± 4.0 and 68.8 ± 3.8 M for the components. The 
final masses strongly depend on both the period and inclination of the binary, which can only 
be measured from the light curve. They derived a period of 3.675 days from the absolute 
values of the radial velocity differences and assumed a circular orbit to derive the masses of 
the components. In their discussion of WR 20a, Rauw et al. (2004) stressed the importance 
and necessity of photometric monitoring of this binary. 

Shara et al. (1991) were the first to make the discovery that WR 20a was a WR star, 
by obtaining spectroscopy of the system. Moffat et al. (1991) obtained UBV photometry, 
however, as Rauw et al. (2004) suggest, it is likely affected by photometric eclipses. WR 20a 
was classified as a candidate binary along with 30 other WR stars in the Vllth Catalogue of 
Galactic Wolf-Rayet Stars (van der Hucht 2001), because of its comparatively weak emission 
lines, which were thought to be diluted by the continuum light of an OB-type companion. 
Rauw et al. (2004) spectroscopically surveyed these candidate WR binaries and WR 20a 
clearly emerged as an interesting object, which they followed up. Due to its favorable location 
in the sky, this object was included in the target list covered by the Optical Gravitational 
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Lensing Experiment (OGLE) photometric survey. We have quickly obtained a precise light 
curve of WR 20a and thus confirm and refine, by measuring the period and inclination, the 
remarkable masses of the WR 20a components. 

In this paper, we present photometry for WR 20a, which clearly shows the presence of 
eclipses. In Section 2 we describe the observations, in Section 3 we present the light curve 
and the analysis and in Section 4 we discuss the results. 

2. OGLE Observations 

WR 20a (= SMSP2 = THA35 - II - 036) is located in the constellation Carina at a = 
10:23:58.0, 5 = -57:45:49 (J2000.0). Observations were carried out with the 1.3 m Warsaw 
telescope at Las Campanas Observatory, Chile, which is operated by the Carnegie Institute 
of Washington. The telescope is equipped with a mosaic CCD camera with 8192x8192 
pixels. WR 20a is very bright, / ~ 11.0, so very short exposures of 10 seconds each were 
used to avoid saturating the CCD. The observations presented here started on May 1, 2004 
and consist of 83 measurements obtained on 17 consecutive nights. 

The absolute calibration of our photometry was obtained on one photometric night 
with a couple of standards from Landolt (1992). The uncertainty in the zero points of the 
photometry is about 0.03 mag. 

The data were reduced using the standard OGLE-III data pipeline, as described by 
Udalski (2003). The photometry was derived using difference imaging analysis. This method 
is the current state-of-the-art for photometric accuracy in crowded fields (Alard 2000). Our 
analysis resulted in a set of time-series photometry in the /-band, presented in Table 1. 
The formal error in the photometry is as low as 0.001 — 0.002 mag, but we have adopted 
a conservative error of 0.01 mag for the differential light curve, y-band images were also 
obtained on one night during maximum. We measured the following magnitudes for WR 20a 
at maximum: V = 13.45 ± 0.025, and / = 10.65 ± 0.05. 

3. Light Curve and Radial Velocity Curve Analysis 

The first goal of our observations was to confirm and refine the orbital period of the 
WR 20a binary. As discussed extensively by Rauw et al. (2004), besides the main period of 
P = 3.675 days, their spectroscopic data also allowed for a much shorter period of 1.293 days, 
which would result in reduced minimum masses of about 25 M for each star in the system. 
A longer period of 4.419 days was also possible, resulting in increased minimum masses of 
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above 80 M for each star. Using the analysis of variance technique of Schwarzenberg-Czerny 
(1989), the most significant period derived from our data is P = 1.843 ± 0.005 days, where 
the error is a conservative 3a estimate. However, we know that the system has to have two 
eclipses, which results in a true orbital period twice as long. We find an ephemeris for the 
primary (deeper) eclipse of 

T prim = 2453124.569 + 3.686 x E, (1) 

thus confirming and refining the orbital period of WR 20a and the remarkable masses of its 
components. 

In addition, we can derive an inclination angle i for the system from our well-sampled 
light curve. To first order, the value of the inclination angle is given by the depth of the 
eclipses (~ 0.4 mag), and we find that the exact value of i only weakly depends on the details 
of the model fit, or even which eclipsing binary model we use to fit. To demonstrate this, 
we first fit the light curve with a simple model of two spherical stars with limb darkening 
(J. Devor 2004, private communication), as well as with the more complex EBOP model 
(Nelson & Davis 1972; Popper & Etzel 1981) and finally, with the Wilson-Devinney (WD) 
code (Wilson & Devinney 1971; Wilson 1979; van Hamme & Wilson 2003) for modeling 
distorted stars. We find that indeed the inclination angle is insensitive to the model, and 
the best fit value with WD is % = 74.5 ± 2.0 deg. 

We ran WD in the overcontact mode (MODE 3), fixing T cfn = 42 000K (since WR 20a 
is intermediate between Mk42 and WR 47c, see Rauw et al. 2004; Crowther & Dessart 1998), 
using linear limb darkening coefficients from van Hamme & Wilson (2003), values for gravity 
darkening exponents and albedos from theoretical values for radiative envelopes, and a mass 
ratio determined from the radial velocity curve (see below). We fit for the inclination i, T c g 2 , 
the luminosity of the primary and the surface potential (Qi = f^) and defined convergence 
to be reached after five consecutive iterations for which the corrections for all adjusted 
parameters were smaller than their respective standard (statistical) errors. The results of 
the fit are given in Table 2. In Figure 1 we show the result of the /—band light curve model 
fit for WR 20a. The uneven eclipse depths suggest slightly different effective temperatures 
and, thus, different spectral types for the components. Using differential photometry we 
can measure a difference between the two stars as small as 1-2%, which is not possible with 
spectral decomposition. However, multi-band photometry is necessary to resolve the effective 
temperature - radius degeneracy. Finally, we derive the following best fit fractional radii 
for both stars: polar radius r polc = 0.34, radius to the inner Lagrangian point r point = 0.40, 
r side = 0.35, and r bac k = 0.37. 
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Armed with the refined period and the exact value of the inclination angle, we re- 
analyzed the radial velocity data of Rauw et al. (2004). Fixing the eccentricity to 0, as the 
light curve confirms, and the period to 3.686 days, we fit the radial velocity data applying 
equal weights and rejecting points with radial velocity measurements smaller than 80 km s -1 . 
These measurements were made by fitting Gaussians to lines that are not resolved and thus 
include errors much larger than 10 — 15 km s _1 and are also subject to non-Keplerian 
effects. Our fit is shown in Figure 2, and the best fit values with their statistical errors 
are presented in Table 3. We note that due to the refined period, as Rauw et al. (2004) 
caution, their two earliest RV epochs are now "flipped". We derive slightly larger velocity 
semi-amplitudes than Rauw et al. (2004), which in turn yield larger masses. Specifically, the 
values are Ki = 362.2 km s _1 and K 2 = 366.4 km s -1 , the mass ratio q = m 2 /mi = 0.99, the 
systemic velocity 7 = 8 km s" 1 , the semi-major axis a sin % = 53 R , and finally the masses 
mi sin 3 i = 74.3 M Q and m 2 sin 3 i = 73.4 M . Knowing the inclination, we calculate the 
masses of the stars to be 83.0 ± 5.0 M and 82.0 ± 5.0 M . These erros include the error in 
the minimum mass and the error in the inclination. 



4. Discussion 

As discussed at length by Rauw et al. (2004), the extreme mass values of WR 20a make 
it a "cornerstone system for future investigations of massive star evolution". We confirm and 
refine these extreme values, by measuring the masses of the components to be 83.0 ± 5.0 M 
and 82.0 ± 5.0 M . If the WN6ha spectral type is confirmed, WR 20a will be one of few 
known WR binaries with a mass ratio near unity. Thus, it becomes very important now to 
obtain high resolution spectra and multiband photometry of WR 20a, corrected for eclipses. 

Such future observations will allow an accurate reddening measurement and further- 
more, a direct, independent and accurate distance determination to the system. This will 
resolve the discrepancy in the distance to Westerlund 2, which ranges from 2.5 — 8 kpc (see 
discussion in Churchwell et al. 2004). Several multiband measurements of WR 20a already 
exist. As measured by us (Section 2), WR 20a has V = 13.45 and / = 10.65 at maximum. 
Our V is different from V = 13.58 measured by Moffat et al. (1991), but their measure- 
ment could have been affected by the eclipses. The Two Micron All Sky Survey (2MASS, 
Skrutskie et al. 1997) has measured J = 8.86, H = 8.08 and K = 7.59 for WR 20a. Recent 
infrared observations of the HII region RCW 49 with the Spitzer Space Telescope (Whitney 
et al. 2004) contain measurements of the binary. Belloni & Mereghetti (1994) have obtained 
ROSAT observations of RCW 49, tentatively identifying WR 20a with one of their point 
sources. However, further observations are needed to determine the extinction to WR 20a. 
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Analogs of WR 20a, possibly with even higher mass components, are bound to exist 
in nearby galaxies, such as M31 and M33. There, they would be first found by variability 
searches such as DIRECT (Stanek et al. 1998; Bonanos et al. 2003) as bright eclipsing 
binaries (about 18th magnitude for My = —6.5), which can then be followed by spectroscopy 
to obtain radial velocities. It is our experience (A.Z. Bonanos et al., in preparation) that 
accurate radial velocities can be obtained for such systems with the existing 6.5-10 meter 
telescopes. The advantage of studying eclipsing binaries in these galaxies would be their well 
known distances and low reddening, as systems with high reddening would be too faint to 
follow-up. 

We conclude that WR 20a truly is a "cornerstone system" that deserves detailed addi- 
tional studies. Not only will studies of the most massive binaries constrain models, but also 
will help understand other problems, such as the nature of mass loss and winds of early type 
stars, as well as the connection of these stars to supernovae and gamma ray bursts (Stanek 
et al. 2003; Maund et al. 2004). WR 20a now overtakes the previous record holders WR 22 
(Rauw et al. 1996; Schweickhardt et al. 1999) and R136-38 (Massey et al. 2002) in the most 
massive star competition, by becoming "the most massive star ever weighed." 

We thank G. Torres for a careful reading of and comments on the manuscript and D. 
Sasselov and J. Devor for useful discussions. We also thank the referee, Eric Gosset, for his 
prompt and careful reading of the manuscript and his useful comments. AZB and KZS were 
partially supported by HST Grant HST-GO-09810.06-A. Support for Proposal number HST- 
GO-09810.06-A was provided by NASA through a grant from the Space Telescope Science 
Institute, which is operated by the Association of Universities for Research in Astronomy, 
Incorporated, under NASA contract NAS5-26555. Partial support to the OGLE project was 
provided with the following grants: Polish KBN grant 2P03D02124 to A. Udalski and NSF 
grant AST-0204908 and NASA grant NAG5-12212 to B. Paczyhski. 
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Fig. 1. — Wilson-Devinney model fit of a near contact binary to the J-band light curve of 
WR 20a. The period is 3.686 days, the eccentricity is and the inclination % = 74.5 deg. 
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Fig. 2. — Radial velocity curve fit to measurements of WR 20a made by Rauw et al. (2004). 
The filled circles correspond to the primary (more massive) star. The new semi-amplitude 
values are Ki = 362 km s" 1 and K 2 = 366 km s -1 and the masses are mx sin 3 i = 74.3 M and 
m 2 sin 3 i = 73.4 M Q , which, combined with i = 74.5 deg yield mass values for the components 
of 83.0 M F) and 82.0 M n . 
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Table 1. /-band Photometry of WR 20a 



HJD 



01 



UT 



May 1, 2004 



May 2, 2004 



May 3, 2004 



May 4, 2004 



May 5, 2004 



May 6, 2004 



May 7, 2004 



May 8, 2004 



May 9, 2004 



May 10, 2004 



2453126.57440 
2453126.57581 
2453126.58207 
2453126.58350 
2453127.57488 
2453127.57883 
2453128.54163 
2453128.54306 
2453128.54659 
2453128.54800 
2453129.58401 
2453129.58545 
2453130.52109 
2453130.52253 
2453131.52484 
2453131.52627 
2453132.55993 
2453132.56136 
2453132.65164 
2453132.65307 
2453132.75066 
2453133.56632 
2453133.56775 
2453133.65146 
2453133.65288 
2453133.72505 
2453133.72649 
2453134.57341 
2453134.57483 
2453135.45134 
2453135.45275 
2453135.45523 
2453135.45666 
2453135.50139 
2453135.50284 
2453135.50518 
2453135.50660 



10. 


.932 


0, 


.01 


10. 


.931 





.01 


10. 


.920 


0. 


.01 


10. 


.925 





.01 


10. 


.657 





.01 


10. 


.658 


0. 


.01 


10. 


.804 


0. 


.01 


10. 


.808 


0. 


.01 


10. 


,811 


0. 


.01 


10. 


.797 


0. 


.01 


10. 


.682 


0, 


.01 


10. 


,692 


0. 


.01 


10. 


.723 


0, 


.01 


10. 


,716 


0. 


.01 


10. 


,711 





.01 


10. 


.706 


0. 


.01 


10. 


.687 


0. 


.01 


10. 


.688 


0. 


.01 


10. 


,672 





.01 


10. 


,677 


0. 


.01 


10. 


.641 


0, 


.01 


10. 


,849 





.01 


1 (1 




u. 


m 


10. 


,938 


0. 


.01 


10. 


,952 





.01 


11. 


,029 


0, 


.01 


11. 


,010 


0, 


.01 


10. 


,655 


0, 


.01 


10. 


,658 





.01 


10. 


,918 


0, 


.01 


10. 


,911 





.01 


10. 


,909 





.01 


10. 


,920 


0, 


.01 


10. 


,963 


0, 


.01 


10. 


,976 





.01 


10. 


,976 
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Table 1 — Continued 



UT HJD I 



May 11, 2004 



May 12, 2004 



May 13, 2004 



May 14, 2004 



May 15, 2004 



2453135.56331 
2453135.56473 
2453135.62564 
2453135.62707 
2453135.62904 
2453135.63045 
2453135.72198 
2453135.72339 
2453136.48681 
2453136.48824 
2453137.46474 
2453137.46618 
2453137.50263 
2453137.50407 
2453137.53984 
2453137.54127 
2453137.58009 
2453137.58149 
2453137.62247 
2453137.62389 
2453138.46632 
2453138.46773 
2453138.58603 
2453138.58744 
2453139.48039 
2453139.48180 
2453139.54939 
2453139.55083 
2453139.60100 
2453139.60241 
2453139.63839 
2453139.63982 
2453140.46213 
2453140.46355 
2453140.53695 
2453140.53839 
2453140.58621 
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Table 1 — Continued 



UT 


HJD 


I 


0i 




2453140.58762 


10.670 


0.01 


May 16, 2004 


2453141.55358 


10.726 


0.01 




2453141.55499 


10.721 


0.01 


May 17, 2004 


2453142.46556 


10.705 


0.01 




2453142.46700 


10.700 


0.01 




2453142.52384 


10.710 


0.01 




2453142.52525 


10.717 


0.01 




2453142.57798 


10.725 


0.01 




2453142.57940 


10.725 


0.01 



Table 2. Light Curve Parameters for WR 20a 



Parameter 




Value 


Period, P 




3.686 ± 0.01 days 


Primary eclipse, T pr i m 




2453124.569 


Inclination, i 




74.5 ± 2.0 deg 


Eccentricity, e 




(fixed) 


Effective Temperature, 


T e ffl 


42000 K (fixed) 


Effective Temperature, 


T e ff2 


40300 ± 1000 K 


Surface Potential (f2i = 


= fi 2 ) 


3.92 ± 0.03 


Radius, r po i 




18.7 ± 0.3 R 


Radius, r poin t 




22.0 ± 0.3 R 


Radius, r sido 




19.3 ± 0.3 R Q 


Radius, r bac k 




20.4 ± 0.3 R© 



Table 3. Radial Velocity Curve Parameters for WR 20a 



Parameter 


Value 


Systemic velocity, 7 


8.0 ± 5.0 kms- 1 


Semi-amplitude, Ki 


362.2 ± 8.0 kms" 1 


Semi-amplitude, K2 


366.4 ± 8.0 kms" 1 


Semi-major axis, a sin i 


53.0 ± 1.0 R Q 


Mass Ratio, q 


0.99 ± 0.03 


Mass, mi sin 3 i 


74.3 ±4.OM 


Mass, m 2 sin 3 i 


73.4 ±4.0M Q 



